Multirate integration of axisymmetric step-flow equations
نویسندگان
چکیده
We present a multirate method that is particularly suited for integrating the systems of Ordinary Differential Equations (ODEs) that arise in step models of surface evolution. The surface of a crystal lattice, that is slightly miscut from a plane of symmetry, consists of a series of terraces separated by steps. Under the assumption of axisymmetry, the step radii satisfy a system of ODEs that reflects the steps’ response to step line tension and step-step interactions. Two main problems arise in the numerical solution of these equations. First, the trajectory of the innermost step can become singular, resulting in a divergent step velocity. Second, when a step bunching instability arises, the motion of steps within a bunch becomes very strongly stable, resulting in “local stiffness”. The multirate method introduced in this paper ensures that small time steps are taken for singular and locally stiff components, while larger time steps are taken for the remaining ones. Special consideration is given to the construction of high order interpolants during run time which ensures fourth order accuracy of scheme for components of the solution sufficiently far away from singular trajectories.
منابع مشابه
Multirate Undrained Splitting for Coupled Flow and Geomechanics in Porous Media
We consider a multirate iterative scheme for the quasi-static Biot equations modelling the coupled flow and geomechanics in a porous medium. The iterative scheme is based on undrained splitting where the flow and mechanics equations are decoupled with the mechanics solve followed by the pressure solve. The multirate scheme proposed here uses different time steps for the two equations, that is, ...
متن کاملAnalytical study of flow field and heat transfer of a non-Newtonian fluid in an axisymmetric channel with a permeable wall
In this study, the momentum and energy equations of laminar flow of a non-Newtonian fluid are solved in an axisymmetric porous channel using the least square and Galerkin methods. The bottom plate is heated by an external hot gas, and a coolant fluid is injected into the channel from the upper plate. The arising nonlinear coupled partial differential equations are reduced to a set of coupled no...
متن کاملAixsymmetric Stagnation Point Flow of a Viscous Fluid on a Moving Cylinder with Time Dependent Axial Velocity
The unsteady viscous flow in the vicinity of an axisymmetric stagnation point of an infinite moving cylinder with time-dependent axial velocity is investigated. The impinging free stream is steady with a strain rate k. An exact solution of the Navier-Stokes equations is derived in this problem. A reduction of these equations is obtained by use of appropriate transformations. The general self-si...
متن کاملRational Chebyshev Collocation approach in the solution of the axisymmetric stagnation flow on a circular cylinder
In this paper, a spectral collocation approach based on the rational Chebyshev functions for solving the axisymmetric stagnation point flow on an infinite stationary circular cylinder is suggested. The Navier-Stokes equations which govern the flow, are changed to a boundary value problem with a semi-infinite domain and a third-order nonlinear ordinary differential equation by applying proper si...
متن کاملA multirate time integrator for regularized Stokeslets
The method of regularized Stokeslets is a numerical approach to approximating solutions of fluid-structure interaction problems in the Stokes regime. Regularized Stokeslets are fundamental solutions to the Stokes equations with a regularized point-force term that are used to represent forces generated by rigid or elastic object interacting with the fluid. Due to the linearity of the Stokes equa...
متن کامل